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Abstract 

We suggest that if the astrophysical site for r-process nucleosynthesis in the early Galaxy is confined to 
a narrow mass range of Type II supernova (SN II) progenitors, with a lower mass limit of = 2OM0, 
a unique feature in the observed distribution of [Ba/Mg] vs. [Mg/H] for extremely metal-poor stars can 
be adequately reproduced. We associate this feature, a bifurcation of the observed elemental ratios into 
two branches in the Mg abundance interval — 2.7 < [Mg/H] < —2.3, with two distinct processes. The first 
branch, which we call the -branch, is associated with the production of Ba and Mg from individual 
massive supernovae. The derived mass of Ba synthesized in SNe II is 8.5 x lO^^M© for Mms = 2OM0, 
and 4.5 x lO-^M© for M^s = 25M0. We conclude that SNe II with M^s ^ 2OM0 are the dominant 
source of r-process nucleosynthesis in the early Galaxy. An SN-induced chemical evolution model with this 
Mms-dependent Ba yield creates the ^/-branch, refiecting the different nucleosynthesis yields of [Ba/Mg] for 
each SN II with ^ 2OM0. The second branch, which we call the "z" -branch, is associated with the 
elemental abundance ratios of stars which were formed in the dense shells of the interstellar medium swept 
up by SNe II with M^s < 2OM0 that do not synthesize r-process elements, and applies to stars with 
observed Mg abundances in the range [Mg/H]< —2.7. The Ba abundances in these stars refiect those of the 
interstellar gas at the (later) time of their formation. The existence of a [Ba/Mg] i-branch strongly suggests 
that SNe II which are associated with stars of progenitor mass Mms ^ 2OM0 are infertile sources for the 
production of r-process elements. We predict the existence of this i-branch for other r-process elements, 
such as europium (Eu), to the extent that their production site is in common with Ba. 
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supernovae: general — supernova remnants 



1. INTRODUCTION 

Truran (1981) suggested that the observed presence of 
elements heavier than iron (Fe) in extremely metal-poor 
stars is due to the nucleosynthesis products of rapid neu- 
tron capture reactions (the r-process). Europium (Eu), 
for which about 97% of its abundance in the solar sys- 
tem is of pure r-process origin (Kappler, Beer, & Wis- 
shak 1989), has therefore been considered as one of the 
most useful elements for locating the astrophysical sites 
of r-process nucleosynthesis. Assuming that such sites 



would most likely be identified with Type II supernovae 
(SNe II), many authors have attempted to constrain the 
applicable mass range of SN II progenitors from an ob- 
served [Eu/Fe] vs. [Fe/H] relation for metal-poor stars. 
However, a satisfactory mass range for the progenitor 
stars has not yet been agreed upon. Mathews, Bazan, & 
Cowan (1992) suggested a mass range of = 7 — 8M0 
for the astrophysical r-process site, while Travaglio et 
al. (1999) supported a somewhat higher mass range, 
= 8 — 1OM0. Most recently, Ishimaru & Wanajo 
(1999) concluded that the [Eu/Fe] versus [Fe/H] relation 
alone is unable to distinguish an r-process site with pro- 
genitor masses Mms = 8 — IOM0 from that resulting from 
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a site with M^s > SOM©. 

The discovery of numerous Galactic halo stars with 
extremely low metal abundances from the ongoing HK 
survey of Beers and colleagues (Beers, Preston, & Sheet- 
man 1992; Beers 1999), in particular those in the abun- 
dance range — 4.0 < [Fe/H] < —2.5, have opened the door 
for much more detailed investigations of the elemental 
production from SNe II in the early Galaxy. Shigeyama 
& Tsujimoto (1998) recently argued that the elemental 
abundance patterns observed in the atmospheres of ex- 
tremely metal-poor stars retain those produced by indi- 
vidual SNe II (see also Audouze & Silk 1995), and thus 
can be used to estimate the heavy element yields of the 
first generation of stars in the Galaxy. They also argued, 
on the basis of theoretical SN progenitor models, that 
the synthesized mass of magnesium, MMg, as a function 
of Mms, is likely to be far more certain than the pre- 
dicted mass of iron produced, Mpe, which is dependent 
on the mass cut chosen in the stellar core. These consid- 
erations indicate that the observed [Eu/Mg] vs. [Mg/H] 
relation in extremely metal-poor stars should give a much 
firmer constraint on the r-process site. According to Tsu- 
jimoto & Shigeyama (1998), stars with the lowest value 
of [Mg/H]= -2.7 (CS 22892-052) that exhibit absorption 
lines of Eu in their spectra in the sample of McWilliam 
et al. (1995) set the lower limit to ^ 2OM0, above 
which SNe II produce the r-process elements. 
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Fig. 1. The correlation of [Ba/Mg] with [Mg/H] for 
metal-poor stars. Filled circles show the data of 
McWiUiam (1998). Two stars (CS 22898-027, CS 
22947—187) in his sample are excluded from the 
plots because their Ba abundances may be contam- 
inated by s-process material accreted from binary 
companions (McWilliam 1998). Other r-process 
data (McWilliam et al. 1995) are also shown by 
different symbols indicated in the figure. 

For most stars of extremely low metallicity ([Fe/H] 
, — 2.0), obtaining an abundance measurement of Eu 



with 2.5-4m-class telescopes is a difficult observational 
task. For example, McWilliam et al. (1995) only ob- 
tained Eu abundances for 11 of the 33 metal-poor stars in 
their sample. This incompleteness implies that one might 
miss valuable insights to r-process production mecha- 
nisms by concentrating exclusively on observations of Eu. 
McWilliam (1998) showed that the [Ba/Eu] ratios for 
stars with —3^ [Fe/H] ^ —2 are consistent with pure r- 
process ratios, and concluded that Ba in extremely metal- 
poor stars is also of r-process origin. Burris et al. (1999) 
arrived at a similar result. 

In Figure 1, which uses the data of McWilliam et 
al. (1995) and McWilliam (1998), the [Ba/Mg] values for 
a sample of metal-poor stars are plotted against [Mg/H], 
together with four heavier r-process elements. La, Ce, 
Nd, and Eu. It is indeed remarkable that there ap- 
pears a nearly vertical boundary at [Mg/H]~ —2.5, with 
[Ba/Mg] spanning the range from [Ba/Mg] = —2 to +0.6, 
whereas a horizontal line of [Ba/Mg] ~ —1.4 emerges from 
[Mg/E]^ -2.5 down to -3.7. The data of Ryan, Norris, 
& Beers (1996) for seven of their extremely metal-poor 
stars show a similar trend in the [Ba/Mg] — [Mg/H] plane. 

In §2 we discuss the r-process sites which might be 
consistent with the vertical [Ba/Mg] — [Mg/H] boundary, 
based on the models proposed by Shigeyama & Tsuji- 
moto (1998) and Tsujimoto & Shigeyama (1998). In 
§3 the entire range of observed [Ba/Mg] — [Mg/H] abun- 
dances in halo stars is used to explore the early stages 
of inhomogeneous chemical evolution in the Galaxy in 
the context of the SN-induced star-formation model de- 
scribed by Tsujimoto, Shigeyama, & Yoshii (1999, here- 
after TSY). 



2. PRODUCTION SITE FOR r-PROCESS EL- 
EMENTS 

Figure 2a illustrates the hypothesis we seek to defend, 
that the metal-poor stars of the Galaxy populate two sep- 
arate branches in the [Ba/Mg] — [Mg/H] plane. The first 
branch extends rightward of the vertical boundary which 
begins at ([Ba/Mg], [Mg/H]) = (-2.0, -2.3) and ends at 
([Ba/Mg], [Mg/H]) = (+0.6,-2.7). The second branch is 
horizontal from [Mg/H]~ —2.5 to —3.7 at a constant 
value of [Ba/Mg]- -1.4. 

Other r-process elements, such as La, Ce, Nd, and 
Eu, shown in Figure 1, populate the first branch iden- 
tified from the behavior of [Ba/Mg] in Figure 2, but we 
presently lack measurements of their abundances in the 
most metal-poor stars, which are required in order to 
confirm the presence of the second branch. Tsujimoto 
& Shigeyama (1998) derived the Eu yield as a function 
of SN II progenitor mass, and showed that average pre- 
dicted value, weighted by the Salpeter initial mass func- 
tion (IMF), successfully reproduces a plateau value of 
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Fig. 2. (a) The correlation of [Ba/Mg] with [Mg/H] 
for metal-poor stars; same as in Fig.l. The thick 
line is the left-side boundary to the y-branch rep- 
resenting the mass range of the r-process site. (6) 
SN II progenitor mass Mms plotted as a function of 
[Mg/H] inside the shell swept up by the SNR. The 
stellar lifetimes are also shown for four reference 
masses of Mms = 12, 20, 25, and SOM©. Dashed 
lines in the lower and upper panels illustrate how 
we assign Mms to each star along the boundary. 



[Eu/Fe] seen at -2 < [Fe/H] < - 1. Thus, hereafter, we 
refer to the first branch as the ^-branch, where the letter 
"i/" stands for its origination in individual SNe II yields. 

As shown in Figure 2a, the ^/-branch is confined to stars 
with [Mg/H]> -2.7; there is no star with [Mg/H] < -2.7 
that belongs to the ^/-branch. On the other hand, 
Shigeyama & Tsujimoto (1998) showed that the metal- 
licity []VIg/H] of stars born from an SN remnant (SNR) is 
well approximated by the average [IVIg/H] inside the shell 
swept up by the SNR. This gives a relation between the 
metallicity []VIg/H] of stars and the mass M^s of SN II 
progenitor as shown in Figure 2b. As a consequence, it is 
indicated that only SNe II with Mms ^ 20Mq produce Ba 
via the r-process. If we assume that the vertical bound- 
ary to the ^/-branch has a one-to-one correspondence to 
the Ba yield from individual SNe II, the progenitor mass 



range is confined to M^s = 20 — 25Mq. It is straightfor- 
ward to derive the Ba yield from the observed [Ba/Mg] - 
value along the vertical boundary and the synthesized Mg 
mass, leading to Msa = 8.5 x IQ-^Mq for M^s = 2OM0, 
and MBa = 4.5 x 10~^Mq for M^s 25M0. 

We now consider the origin of the second branch of 
[Ba/Mg], in the range of -3.7 < [Mg/H] < - 2.7. This 
[Mg/H] range corresponds to M^s = 12 - 20 Mq, for 
which SNe II, according to our present models, do not 
produce significant amounts of Ba. Therefore, our hy- 
pothesis is that the Ba abundances which are observed 
in the atmospheres of these stars come only from the 
interstellar matter (ISM) that was enriched in Ba by 
the preceding SNe II (with M^s = 20 - 25M0), and 
that was swept up in the shells by later SNe II with 
Mms = 12 — 2OM0. Thus, hereafter we refer to this 
branch, as the i-branch, where the letter "i" stands for 
an ISM origin. We note that the [Ba/Mg] value in 
this branch should always be below its plateau value 
[Ba/Mg] ~ —0.6 at higher metallicities. We also predict 
from the above argument that stars born from the shell 
swept up by SNe II with Mms > 25M0 would form an- 
other i-branch at [Mg/H] ^ — 2.5, evidence for which is 
not seen in the sample of McWilliam et al. (1995). 



3. CHEMICAL EVOLUTION OF i^-PROCESS 
ELEMENTS IN THE GALACTIC HALO 

In this section, we discuss the chemical evolution of Ba 
and Eu in the Galactic halo, based on the formulation 
presented in TSY. The essence of the TSY model is that 
the star-forming process is assumed to be confined to 
separate clouds which make up the entire halo, and that 
the chemical evolution in these clouds proceeds through 
a successive sequence of SN II explosion, shell formation, 
and resulting star formation. Details can be found in 
TSY. 

An overview of our model applied to the Galactic halo 
is as follows: - (1) The metal- free Population HI stars 
(Pop HI) form by some (as yet unspecified) mechanism 
in primordial-composition gas clouds of the Galactic halo. 

(2) The most massive stars among them explode as Pop 
HI SNe II, which trigger a series of star formation events. 

(3) Star formation terminates when SNRs become unable 
to the formation of dense shells. (4) Roughly 90 % of the 
cloud mass remains unused in star formation, and may 
fall onto the still-forming Galactic disk. 

The free parameters in our model are the mass fraction 
xui of metal-free Pop HI stars initially formed in each 
cloud, and the mass fraction e of stars formed in a dense 
shell swept up by each SNR. A value of xm sensitively de- 
termines the level of [Ba/Mg] for the i-branch, so we take 
xiii = 2.5 X 10~^ to be consistent with the observed value 
[Ba/Mg]=-1.4at [Mg/H]- -3.7 (McWilliam 1998). We 
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take e = 4.3 X 10 ^ to reproduce the observed [Fe/H] 
distribution function of halo field stars for [Fe/H]< — 1 
(TSY). 
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Fig. 3. (a) The color-coded frequency distribution of 
the long-hved stars in the [Ba/Mg] — [Mg/H] plane, 
convolved with a Gaussian having (7=0.2dex for 
[Ba/Mg] and cr=0.1dex for [Mg/H]. The symbols 
represent the data taken from McWilliam 1998 
(filled squares), Ryan et al. 1996 (open circles), 
Burris et al. 1999 (filled triangles), Luck Sz Bond 
1985 (filled circles). Luck & Bond 1981 (open trian- 
gles), Magain 1989 (open squares), and Steenbock 
1983 (crosses). For comparison, the result of the 
conventional one-zone chemical evolution model is 
shown by dashed line assuming the r-process site 
of 8 — 10 Mq supernovae, which is different from 
our proposed site of 20 — 25 Mq supernovae. (b) 
The same as (a) but for Eu. The dashed curve en- 
closes the predicted region in which i-branch stars 
must exist. 

Adopting the result of §2, that SNe II with = 
20 — 25M0 are the dominant site for r-process nucle- 
osynthesis, and using the derived Ba yield dependent on 
Mms, we calculate the Ba evolution at early epochs of the 
Galactic halo. For later evolution, which pertains to the 
abundance range — 2 ;S []VIg/H] ^ —1, we add the contri- 
bution from 5-process nucleosynthesis, assuming that the 
production site for early-epoch s-process elements is 2 — 3 



Mq AGB stars (Burris et al. 1999). The stellar IMF is 
assumed to be of the Salpeter form. The upper and lower 
mass limits for stars that explode as SNe II are taken to 
be 5OM0 and I2M0, respectively. 

Figure 3a is a color-coded predicted frequency distri- 
bution of stars in the [Ba/lVIg] — []VIg/H] plane, normal- 
ized to unity when integrated over the entire area (see 
the color bar for the scale). In order to enable a di- 
rect comparison with the data, the frequency distribu- 
tion has been convolved with a Gaussian with a = 0.2 
dex for [Ba/Mg] and a = 0.1 dex for [Mg/H]. As already 
pointed by Burris et al. (1999), there is a systematic dif- 
ference in the derived Ba abundance between Burris et 
al. (1999) and McWilliam (1998). Burris et al. (1999) 
obtain [Ba/H] ratios which are higher, in the mean, by 
0.46 dex for the seven stars in common with the data of 
McWilham (1998). We therefore decrease their [Ba/H] 
ratios by 0.46 dex to enable a consistent comparison with 
the data of McWilliam (1998) in Figure 3a. The [Ba/Mg] 
ratio "plateau" in the range of — 2 ;S [Mg/H] ^ — 1 is ob- 
viously higher than that of z-branch ([Mg/H] ^ —2.5), 
which implies a significant contribution from 5-process 
nucleosynthesis. For comparison, we show by the dashed 
line in this figure the conventional one-zone chemical evo- 
lution model assuming that the production sites for r- 
process and s-process are 8 — 10 Mq supernovae and 2 — 3 
Mq AGB stars, respectively. It is clear that the conven- 
tional model predicts a much smaller [Ba/Mg] than is 
consistent with the observed z-branch at [Mg/H] ^ —2.7. 

To summarize our hypothesis: The earliest stars, stars 
that were formed from the shells swept up by SNe II 
with = 20 — 25M0, first create the vertical bound- 
ary seen in Figure 2a. Thereafter, the ^/-branch devel- 
ops at higher [Mg/H] and converges to the plateau level 
of [Ba/Mg] at higher metallicities, as described in de- 
tail in TSY. At the same time, the subsequently formed 
stars, arising from the SNRs with M^s = 12 — 2OM0, do 
not produce Ba, and therefore have the lowest value of 
[Ba/Mg] at [Mg/H] < - 2.5. These stars populate the 
lower envelope of the z-branch, then fill in this branch 
at lower [Mg/H], but gradually increasing [Ba/Mg]. As 
the Mg abundance in the ISM becomes dominant over 
that ejected from SNe II, the z-branch develops at higher 
abundances, increasing both [Mg/H] and [Ba/Mg]. 

Because the upper mass limit in our present model for 
the Ba production site, 25M0, is set below the canon- 
ical upper limit ~ 5OM0 for stable stellar masses, we 
predict the existence of another z-branch in the range of 
—2.3 ^ [Mg/H] ^ — 1.5, which is made up of stars aris- 
ing from the SNRs with M^s = 25 - 5OM0. While this 
branch does not show up in the sample of McWilliam 
(1998), the data from Luck & Bond (1981,1985), inter- 
estingly, support such a prediction (see filled and open 
circles in Fig. 3a). If we increase the upper mass limit to 
5OM0 for the Ba production site, a small amount of Ba is 
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produced in SNe with M^s ^ 25M0, which only slightly 
contaminates the z-branch stars. The results of our calcu- 
lations with an upper mass limit of Mms = 25M0 remain 
almost unchanged, because of the assumed declining IMF 
at high masses. 

The value of xm = 2.5 x 10~^ adopted here can be 
converted into a probability of finding Pop III stars, 
Pill = 2.5 X 10~^, defined as the expected number of 
Pop III stars divided by the total number of the long- 
lived stars in a sample (see TSY). However, we regard 
this xiii-value as an upper bound, because the predicted 
[Ba/Mg] vs. [Mg/H] distribution with smaller xm pro- 
vides an adequate fit to the observed distribution. A 
much more strict constraint on pm could be obtained if 
the sample of extremely metal-poor stars with measured 
Ba abundances is increased by a factor of ten. 

As a final exercise, we also calculate the Eu evolu- 
tion in the early Galaxy using the same models and as- 
sumptions we have employed for the prediction of Ba 
evolution. Figure 3b shows the color-coded frequency 
distribution of stars in the [Eu/Mg] — [Mg/H] plane, after 
convolving with a = 0.2 dex for [Eu/Mg] and a = 0.1 
dex for [Mg/H]. The Eu production site must be the 
same as for Ba, or at least share a part of SNe II sites 
where Ba is produced, because the rapid n-capture pro- 
cess cannot synthesize Eu without producing Ba. There- 
fore, we assume that Eu is also produced by SNe II 
with Mms = 20 - 25M0. We scale the Ba yield by 
the pure r-process value [Ba/Eu]=— 0.72 (Wisshak, Voss, 
& Kapeler 1996) and obtain a predicted Eu yield of 
Meu = 1.3 X lO-^M© for M^s = 20Mq, and Meu = 
7.0 X 1O"^M0 for Mms = 25M0. Our calculation pre- 
dicts the existence of a Eu z-branch at [Mg/H] ^ — 2.5, 
which must be confirmed by future measurements of Eu 
abundances for extremely metal-poor stars. □ 



4. CONCLUSION 

We have suggested the existence of two distinct corre- 
lations of [Ba/Mg] with [Mg/H] for metal-poor stars in 
the Galaxy, the y- and z-branches, separated at [Mg/H] 
~ —2.7. These branches cross one another almost per- 
pendicularly, and form an arrow-like frequency distribu- 
tion of stars in the [Ba/Mg] — [Mg/H] plane. The verti- 
cal boundary to the ^/-branch extending from [Ba/Mg] = 
—2.0 to +0.6 refiects the different nucleosynthesis con- 
tributions of [Ba/Mg] for each SN II of progenitor mass 
in a narrow range, M^s = 20 — 25M0. The horizontal 
z-branch extending from [Mg/H]= —2.7 to —3.7 is pop- 
lin TSY, we used the data of Ryan et al. (1996) for determining 
the Eu yield. These observations exhibit a much lower [Eu/Mg] 
ratio than that of McWilliam et al. (1995). As a consequence, the 
entire population of SNe II with Mms = 12 — 50Mq may be the Eu 
production site. 



ulated by stars born from the shells swept up by SNe H 
with Mms = 12 — 2OM0 that do not synthesize the r- 
process elements. Therefore, the Ba abundances of the 
z-branch stars reflect those in the ISM at the time when 
such stars form. 

The i-branch provides strong evidence for the existence 
of a lower mass limit of the SN progenitors for the site 
of r-process nucleosynthesis. The upper mass limit for 
the r-process site cannot be determined well from the 
current observations. Because of the declining IMF, SNe 
II with Mms ^ 25M0 contributes little to the total Ba 
production in the early Galaxy. In other words, in our 
formulation, the r-process site is dominated by SNe II 
with Mms ~ 20 - 25M0. 

The z-branch should not exist for a- and Fe-peak 
heavy elements formed in the early Galaxy, because 
the entire range of SNe II progenitor masses will con- 
tribute to substantial amounts of such elements. How- 
ever, the z-branch is expected to exist, though is not 
confirmed as of yet, in the case of the r-process ele- 
ments such as La, Ce, Nd, and Eu. It is therefore im- 
portant to measure the abundances of such elements for 
the stars CS 22878-101, CS 22897-008, CS 22891-200, 
CS 22891-209, CS 22948-066, CS 22952-015, and CS 
22950-046 in the sample of McWilliam et al. (1995), 
which actually form the Ba z-branch. 
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ture, and Sports in Japan. 



References 

Audouze, J., & Silk, J. 1995, ApJ, 451, L49 
Beers, T. C, in The Third Stromlo Symposium: The 
Galactic Halo, eds. B. Gibson, T. Axelrod, & M. Put- 
man (ASP: San Francisco), 165, 206 
Beers, T. C, Preston, G. W., & Shectman, S. A. 1992, 
AJ, 103, 1987 

Burris, D. L., Pilachowski, C. A., Armandroff, T. E., 
Sneden, C, Cowan, J. J., & Roe, H. 1999, ApJ, in 
press 

Ishimaru, Y., & Wanajo, S. 1999, ApJ, 511, L33 
Kappler, F., Beer, H., & Wisshak, K. 1989, 

Rep.Prog.Phys., 52, 945 
Luck, R. E., & Bond, H. E. 1981, ApJ, 244, 919 
Luck, R. E., & Bond, H. E. 1985, ApJ, 292, 559 
Magain, P. 1989, A&A, 209, 211 

Mathews, G. J., Bazan, G., & Cowan, J. J. 1992, ApJ, 
391, 719 



6 



T. Tsujimoto et al. 



McWilliam, A., Preston, G. W., Sneden, C, & Searle, L. 

1995, AJ, 109, 2757 
McWilliam, A. 1998, AJ, 115, 1640 
Ryan, S. G., Norris, J. E., & Beers, T. C. 1996, ApJ, 471, 

254 

Shigeyama, T., & Tsujimoto, T. 1998, ApJ, 507, L135 

Steenbock, W. 1983, A&A, 126, 325 

Travaglio, C., Galli, D., Gallino, R., Busso, M., Ferrini, 

F., & Straniero, O. 1999, ApJ, 521, 691 
Tsujimoto, T., & Shigeyama, T. 1998, ApJ, 508, L151 
Tsujimoto, T., Shigeyama, T., & Yoshii, Y. 1999, ApJ, 

519, L64 

Wisshak, K., Voss, F., & Kapeler, F. 1996 in Proc. 8th 
Workshop on Nuclear Astrophysics, ed. W. Hillebrandt 
& E. Miiller (Garching: MPIA), 16 



